Heterogeneous diffusion in a reversible gel 
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We introduce a microscopically realistic model of a physical gel and use computer simulations to 
study its static and dynamic properties at thermal equilibrium. The phase diagram comprises a 
sol phase, a coexistence region ending at a critical point, a gelation line determined by geometric 
percolation, and an equilibrium gel phase unrelated to phase separation. The global structure of the 
gel is homogeneous, but the stress is only supported by a fractal network. The gel dynamics is highly 
heterogeneous and we propose a theoretical model to quantitatively describe dynamic heterogeneity 
in gels. We elucidate several differences between the dynamics of gels and that of glass-formers. 
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Although gels are commonly used in everyday life 
they continue to offer fundamental challenges to research. 
Their physics is determined by a wide window of length- 
scales, from the molecular size of particles in the solvent 
to macroscopic structures, and by a similarly broad range 
of timescales: Gels are "complex" fluids Of partic- 
ular interest are physical gels which are typically made 
of molecules forming a stress-sustaining network, with 
links that have a finite lifetime, as opposed to chemical 
gels where junctions are permanent and properties follow 
directly from geometry. The transient character of the 
network in physical gels results in a complex interplay 
between structure and dynamics, leading to non-trivial 
flow properties. Here we propose a model of a reversible 
physical gel which is microscopically realistic (we are in 
fact inspired by one particular material), and specifically 
design a hybrid Monte Carlo / molecular dynamics nu- 
merical approach to successfully bridge the gap between 
microscopic details and macroscopic observations, while 
offering deep insight on the nature of physical gels. 

Inspired by recent experimental work on gelation, a 
variety of "minimal" models have recently been studied 
to elucidate the connection between gelation and seem- 
ingly related phenomena: Geometric percolation 0, 0|, 
glass transition 0, [H, 0] , phase separation 0, @] . Detailed 
experiments performed with colloidal particles with tun- 
able interactions Q revealed that a non-trivial interplay 
between phase separation and kinetic arrest may produce 
gel-like structures. Associating polymers constitute an- 
other well-studied example of reversible gels [l| . In that 
case, gels can be obtained far from phase separation, pro- 
ducing viscoelastic materials with highly non-linear rhe- 
ological properties that are not well understood 1(| 11 1 . 
In many cases, a close similarity between gelation and 
glass formation is reported [l|. We explain below this 
similarity but discuss also important differences. 

Our model is inspired by a material described in 
Ref. ll|. It is a microemulsion of stable and monodis- 



pcrsc oil droplets in water mixed with tclcchclic poly- 
mers, i.e. long hydrophilic chains with hydrophobic end 
caps. A polymer can form a loop around a single droplet, 
or, more interestingly, a bridge between two droplets. 
Figure [1] is a snapshot taken from our simulations show- 
ing droplets and bridging polymers. For sufficiently high 
polymer concentrations, a percolating network can be 
formed (shown in light gray) and the system becomes a 
soft solid. However, thermally activated extraction of the 
hydrophobic heads leads to a slow reorganization of the 
network structure, and the material eventually flows at 
long times 12| . This material is interesting because func- 
tionality, lifetime of the bonds, volume fraction, strength 
of the networks can all be adjusted independently, which 
is not always possible in attractive colloids @ , or in pre- 
vious model systems, unless specific adhoc assumptions 
are made [3, IS IS • Modeling such a complicated self- 
assembly is a challenge because of the wide range of scales 
involved. In our model we neglect the solvent and include 
polymers and droplets as the elementary objects. More- 
over, since the internal dynamics of the polymers is much 
faster than the gel dynamics, we coarse-grain the poly- 
mer description and only retain their effect as links in- 




FIG. 1: (a) Structure of the physical gel for cj> = 0.2, R = 2 
and TV = 10 4 . A percolating (light gray) cluster of droplets 
connected by telechelic (red) polymers, through which the 
remaining (dark blue) droplets can diffuse, see zoom in (b). 
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ducing an effective entropic interaction between the two 
droplets they connect. Coarse-graining is a crucial step 
for efficient large scale simulations, not used in previous 
models [l3j |. 

We consider an assembly of N droplets of diameter 
a and mass m, interacting, in the absence of polymers, 
with a pair potential typical of soft spheres, Vi(ry) = 
e i ( cr / r ij) i where is the distance between droplets i 
and j, and e\ an energy scale. The potential is cut off 
and regularized at a finite distance, 2.5 a. In addition, 
N p polymers of maximal extension t can form bridges be- 
tween droplets, or loops. Polymer loops have an energy 
cost eo, but no effect on droplet dynamics. On the other 
hand, bridging polymers induce an entropic attraction 
between connected droplets, which we model using the 
classic FENE form, V 2 (n 3 ) = -e 2 ln [l - (r l0 - a) 2 /I 2 ], 
so that polymers act as springs at small elongation, but 
cannot become longer than I. A configuration is speci- 
fied by the droplets positions and velocities, {ri(t), Vj(t)}, 
and by the polymer N x N connectivity matrix, {CV, }, 
where Cy is the number of polymers connecting droplets 
i and j. Summarizing, the Hamiltonian is thus 



N 

n = J2 (y v?+c iieo +E Wite) + c vMra)] ) • (i) 
j>i 



Simulation proceeds by solving Newton's equations for 
the droplets. Lengthscales are given in units of a, energy 
in units of ei, and times in units of \Jmo 2 We use the 
velocity Verlet algorithm with discretization h — 0.005. 
Simultaneously, we use Monte Carlo dynamics to evolve 
polymers. In an elementary move, a polymer is chosen at 
random, and one of its end caps is moved to a randomly 
chosen neighboring droplet. This proposed move is ac- 
cepted with rate r^minfl, exp(— AV^/T)], where AV2 
is the potential energy change during the move, T is the 
temperature, and ri; n k controls the timescale for polymer 
rearrangements. In experiments Tij n k has an Arrhcnius 
behavior associated to the excitation cost for polymer 
extraction. We set £ — 3.5 a ll|, T = 1, eo = 1, and 
e 2 = 50. We found little influence of e 2 on the phase di- 
agram. The relevant control parameters are the droplet 
volume fraction, <j> — ttN/(6V), where V is the volume, 
the number of polymer heads per droplet, R = 2N p /N, 
and Tii n k, which has no influence on static properties. We 
performed simulations for a wide range of parameters, 
R G [0,18], </> G [0.01,0.3], r Unk G [1,10 4 ], N G [10 3 ,10 4 ]. 

The phase diagram, as obtained after a systematic ex- 
ploration of the control parameter space, is shown in 
Fig. O Its topology is in good_ agreement with exper- 
iments in the range studied 




FIG. 2: Phase diagram for a wide range of volume fraction, cj>, 
and number of polymer heads per droplet, R. Symbols are the 
investigated state points in the sol (O); S e l (O)i an< l phase 
separated (o) regions, the yellow point is the approximate 
location of the critical point. Transition lines are sketched. 



phase [14(. We detect phase coexistence from direct vi- 
sualisation and by measuring the static structure factor, 
S(q) = (N^ 1 J^jk exp^Q ' ( r j — r k)}), which exhibits the 
typical q~ 4 behavior at small q. We observe both nucle- 
ation or spinodal regimes depending on the quench depth. 
Interestingly we find that the kinetics of the phase sepa- 
ration is extremely slow, and that the obtained patterns 
are very similar to the ones observed experimentally in 
short-range attractive colloidal suspensions [§]. 

For (j> > 0.05, a broad gel region exists at thermal equi- 
librium between the sol phase at low-i? and phase sepa- 
ration at large- R, see Fig. [2j In the gel phase, a system- 
spanning cluster of polymer-connected droplets emerges, 
which endows the fluid with viscoelastic properties, see 
Fig. Q] The sol-gel transition coincides with geometric 
percolation of polymer-connected clusters. Our gels are 
homogeneous, i.e. S(q) remains typical of a simple fluid, 
as seen in experiments [ll| . However, the spanning clus- 
ter is highly fractal near percolation, and becomes thicker 
deeper in the gel phase. At percolation we find a distri- 
bution of cluster sizes P(n) ~ n -7 , with 7 w 2.2 com- 
patible with random bond percolation, as seen in other 
systems @, 0, & • Finally the structure of the system 
becomes nontrivial when approaching the critical point, 
where S(q) develops a power law behavior with an expo- 
nent close to -1.5 at small q for <f> c m 0.05, R c w 3.5. 

We now show that the gel phase indeed behaves dy- 
namically as a soft viscoelastic fluid. We have investi- 
gated the dynamics by measuring the self-intermediate 
scattering function, F s (q,t) = (N^ 1 J2j exp[iq • (rj(t) — 
r j(0))])i an d the mean-squared displacement, A 2 (t) 



region, the system resembles a dilute assembly of soft 
spheres, and has simple-liquid properties: This is the sol 
phase. Increasing R increases the effective attraction be- 
tween droplets, so that phase separation occurs at large 
R between a low-^, low-i? phase and a large-0, large- R 



111. In the low-0, low-i? (A^Ej l r jW - r j(0)| 2 )- Figure Ha shows the evolu- 
tion of F s (q,t) from the sol to the gel phase. While re- 
laxation is fast and exponential in the sol phase, a slow 
secondary relaxation emerges at percolation. The final 
decay time varies little in the gel phase, but the height of 
the plateau at intermediate times evolves dramatically. 
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FIG. 3: Self-intermediate scattering function for q = 0.46 (a, 
b) and mean-squared displacement (c, d) for N = 10 3 , (f> = 
0.2. (a, c) present the dynamics for rii n k = 10 3 and several 
values of R across percolation (R p ~ 1.85). Viscoelasticity 
continuously emerges at percolation, (b, d) are for R = 4 and 
different values of THnk> which directly controls the long-time 
decay of F s (q,t), while A 2 (t) remains essentially unchanged. 



A similar behavior is found for the coherent scattering 
function, as in experiments [111 ]. Physically the plateau 
reflects the thermal vibrations of an elastic solid on inter- 
mediate timescale, while long-time decay reflects the flow 
of the system: The system is viscoelastic. In Fig. [3}b we 
show that viscous flow is mostly controlled by rii n k, the 
rate for polymer extraction 0, EH . Flow in this system 
occurs when the percola ting network slowly rearranges 
through polymer moves [I2j. Therefore, gelation corre- 
sponds to the continuous emergence, for increasing poly- 
mer concentration, of a plateau in dynamic functions, 
with an almost constant relaxation timescale, controlled 
by the polymer dynamics. Gelation is thus qualitatively 
different from a glass transition where the plateau height 
remains constant with a dramatic increase of relaxation 
timescales [l6| . Coincidence of gelation and percolation, 
put forward in Q or dispelled in happens whenever 
long-lived bonds make cluster restructuration very slow, 
but does not occur in systems where the bond lifetime is 
short at percolation 

Surprisingly, the mean-squared displacements shown in 
Fig.[3]-c-d appear as poor indicators of the dynamics. The 
comparison between Figs. [3}b-d is in fact quite striking. 
While the final relaxation timescale, r, in F s (q,t) scales 
roughly as riink) the diffusivity, D s , extracted from the 
long-time behavior of A 2 (t) ~ 6D s t is almost constant. 
This is reminiscent of the "decoupling" phenomenon, or 
"breakdown" of the Stokes-Einstein relation, reported in 
supercooled fluids While "fractional" breakdown is 
reported in liquids, D s ~ t~**, with £ in the range 0.7- 
0.9 instead of the normal value f = 1 17] . 
£ « 0, quite an extreme case of decoupling 



In supercooled fluids, decoupling phenomena are com- 
monly attributed to the existence of dynamic heterogene- 
ity, that is, the existence of non-trivial spatio-temporal 
distributions of mobilities. The analogy is confirmed in 
Fig. U where we show distributions of droplet displace- 
ments, G s {r,t) = {N' 1 £\ <f(|r-r ( (t) +r, ; (0)|)). Clearly, 
G s exhibits a bimodal character suggesting coexistence of 
slow arrested droplets and fast diffusing droplets. Qual- 
itatively similar distributions were reported in gels [1, 0] 



and glasses [UJ, |19| . Here, the snapshot in Fig. [T] sug- 
gests an obvious explanation for dynamic heterogeneity. 
At any given time, a system-spanning cluster of droplets 
which behaves as a solid on timescales smaller than 
"Hink coexist with droplets which can more freely diffuse 
through this arrested structure. We quantitatively con- 
firm this interpretation in Fig. 3] where G s is decomposed 
over two families of droplets: G s — caGa + (1 — ca)Gm, 
where A (M) stands for droplets that are arrested (mo- 
bile) at time t = 0, ca being the fraction of droplets 
belonging to the percolating cluster. While the central 
peak in G s is dominated by Ga, the large "non-Gaussian" 
tails are dominated by Gm- 

We now propose an analytic model to describe the dy- 
namic heterogeneity of the gel, which incorporates the 
physical idea of a coexistence of a slow, percolating clus- 
ter of connected droplets and fast, more freely diffusing 
droplets, with a dynamic exchange between the two fam- 
ilies set by polymer moves. Similar physical ideas were 
qualitatively discussed earlier 0, [l8|, [2fl| , but were how- 
ever not exploited within a quantitative model. We define 
g a (r,t), the probability that a droplet makes a displace- 
ment r in a time t provided it belongs to family a during 
the whole time interval [0,t], a = A, M. Let p a (t) be 
the probability that a droplet in family a switches for the 
first time to the complementary family, a, at time t, and 
define P a (t) = J t °° p a {t')dt'. Then we have 

G a (r,t) = P a (t)g a (r,t)+ f dtf [A Q (r, t') o G a (r, t — t')] , 

Jo 

(2) 

where A a (r,t) = p a (t)g a (r,t), and o stands for spa- 
tial convolution. The first term describes droplets which 
persist in the same family between and t, while the 
second term captures family exchanges. To keep the 
model simple, we assume the exchange dynamics to be 
homogeneous, with constant transition rates t~ , so 
p a (t) = exp(— t/T a )/r a . In addition, stationarity implies 
that Ca/ta = (1 — ca)/tm- The model can be solved 
analytically in the Fourier-Laplace domain: 



G„(q, s) 



r a A a (q, s) +-r a A a (q, s)A a (q, s) 
1 - A Q (q, s)A a (q, s) 



(3) 



we find here 
Decoupling 



in gels was reported in different systems [18j. 



To obtain analytic fits to the data we assume a 
Gaussian propagator for mobile droplets, Qm{ v i t) = 
(4nD M t)~ 3 / 2 exp[-r 2 /(4D M t)], where D M is an effec- 
tive diffusivity. We treat the droplets attached to the 
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FIG. 4: Distribution of droplet displacements for <j> = 0.2, 
R = 4, Tii n k = 100 and different times for all droplets (O), 
and its decomposition over droplets that are initially free (o) 
or arrested (□). Lines are from our model, Eq. (0. 

cluster as localized in a bounded region of space of linear 
size a, which reflects the thermal vibrations of the elastic 
solid, ffA(r) = (7ra 2 ) -3 / 2 exp[-r 2 /a 2 ]. From Eq. (J3J) we 
obtain a simple (but lengthy) analytic expression for the 
self-intermediate scattering function F s (q, t) by inverse 
Laplace transform. We numerically invert the Fourier 
transform to get the displacement distributions. The free 
parameters of the model are {ca,Dm, a, Ta}, but the first 
three parameters can be fixed by numerical observations. 
The concentration of droplets connected to the percolat- 
ing cluster, ca ~ 0.94, is directly measured. The value 
Dm ~ 0.3 is evaluated from an intermediate-time fit of 
A 2 f (i) restricted to initially mobile droplets. The local- 
ization length, a ks 1.79, is similarly estimated from the 
plateau in A A (t) restricted to droplets initially belonging 
to the percolating cluster. To get the excellent fits shown 
in Fig.[4j we adjust ta = 2- 10 4 , which coincides well with 
the timescale at which A^(t) becomes diffusive, a physi- 
cally sound definition for the average exchange time. 

We were able to fit data for a wide range of parameters 
and find that our model performs similarly well for other 
state points. From Eq. ©, it is easy to predict that r oc 
ta, while D s = (1 — ca)Dm + cacl 2 /(Ata), in quantitative 
agreement with the decoupling data reported in Fig. [3l 
leading to £ = for large rii n k. The diffusion constant 
is in fact entirely dominated by those droplets which do 
not contribute to viscoelasticity, and is therefore a poor 
indicator of the gel dynamics. These results show that 
dynamic heterogeneity in gels can be stronger than in 
supercooled fluids, but its origin is also much simpler: 
The system structure is heterogeneous, Fig. [TJ while no 
such static structure exists in glasses. 

Although motivated by a specific material, the new 
model for reversible gelation proposed in this work sheds 
light on the microscopic aspects of gelation and the het- 
erogeneous dynamics of gel-forming systems. Moreover, 



our numerical findings motivated a simple yet accurate 
analytic modeling of dynamic heterogeneity, which is 
generally applicable to gels. Although slow and het- 
erogeneous dynamics are superficially reminiscent of the 
physics of supercooled fluids, we discussed several quali- 
tative differences between gels and glasses. 
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for discussions, DYGLAGEMEM, MEyC FIS2005-00791 
and Universidad de Granada for financial support. 
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